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Evaluation of the Energy Distribution Function in Physical Adsorption as 
the Problem of Laplace Transform. II. Spatial Distribution of Adsorption 
Energy Evaluated from the Third Gas—Solid Virial Coefficient 


Wyznaczanie różniczkowego rozkładu energii adsorpcji jako problem transformaty 
Laplace'a. II. Określenie przestrzennego rozkładu energii adsorpcji z trzeciego 
współczynnika wirialnego gaz—ciało stałe 


OnpeneJieHne AucbcpepeHHMaJlbHOTo pacripegeenuus SHeprMM a4ĄCOPÓNMM Kak HpOOJeMBI 
TpaHccpopmanuu JIlanzaca. II. Orrpeqejlenne rpocTpaHcTBeHHoro pacnpeqeJleHuA 
3HepruM aĄ4cOopóHuM M3 TpeTbero BMpHMaJbHOro Ko»dxbuuueuTa ra3—TBepqoe Teo 


INTRODUCTION 


As in the first part of this publication series, we shall consider here 
the simple case, when the state of adsorbate molecule is fully described by 
the coordinates of the center of mass of adsorbate molecule. The evaluation 
of the energy distribution function provides, though quantitative, only 
a total description of gas—solid interactions, without giving any informa- 
tion about the macroscopic structure of adsorbent. Such a knowledge of 
gas—solid interactions is sufficient in the adsorption region, where inter- 
action between adsorbed molecules can be negated. At higher adsorbate 
pressures, however, the contributions from the states in which two, three, 
or more adsorbate molecules are near each other, to the thermodynamics 
properties of adsorption systems, become considerable. Therefore, the full 
knowledge of the gas—solid interactions requires not only the energy 
distribution, but also the spatial distribution of adsorption energies. 

In terms of localized adsorption this means the knowledge of the spatial 
distribution for pairs, triples, etc., of adsorption sites. These spatial distri- 
butions are included in appropriate gas—solid virial coefficients [1]. B, 


" 
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involves the function "os Tj), corresponding to the potential of the force 


between two adsorbate molecules, one located at m, and the other at Tą. 
This is a more complicated situation. Not only the energy distribution 
of sites but the spatial distribution for a pair of sites are required. 

Barker and Everett [2] were the first to investigate the hete- 
rogeneity effects, using the virial description formalism. This problem has 
been considered in the last years by Steele [3] and Hill [4]. 

The purpose of this paper is to propose some methods for investigating 
the spatial energy distribution, on the basis of the real energy distribution 
X(e) evaluated in the first part of this publication series. 


THEORETICAL 


Let us consider the expression for the third gas—solid virial coefficient 
Bss (1, 5] 
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and w."=,"(r12) is the pair interaction energy in bulk phase, Ty is the in- 
termolecular distance. The eq. (1) can be rewritten to the form [6] 
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Above V and V, are the volumes of adsorbate and adsorbed phase, re- 
spectively, defined in the first part of this paper. The eq. (3) follows from 
the fact that out the surface phase Vs, the integrand in the integral (1) 
is equal to zero. Next, the contribution from the states in which the 
functions (1+f;) are near maximum to the integral (1), are ominant, 
when considering the positions within the surface phase V;. Thus, the 
first term (1--f;) (1+f:)ð is about 101° times higher than the second term 
for typical adsorption systems at room temperatures [7]. 
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There have been made many attempts to evaluate the integral (3). 
Apart from difficilties involved by heterogeneity, it is a very complicated 
problem. It is due to the fact that the interactions between two adsorbed 
molecules are perturbed by the presence of a solid surface. 

Steele and Halsey [8] and next Freeman and Halsey 
[9] were the first to evaluate theoretically B;s, but without taking into 
account the third-order interactions. Freeman [10] has pointed out 
experimentally the existence of the third-order effects in the adsorption 
of rare gases on graphitized black carbon. The first quantum-mechanical 
caleulations of the third-order interactions are due to Sinanoglu and 
Prvyver? pt): 

Recent results are due to Wolfe and Sams [12] who have found 
that surface area calculations are rather insensitive to these three- 
-body effects. The problem of evaluating the surface area in the presence 
of these three-body effects has been recently considered by Rudziń- 
ski et al. [13] by using the more realistic three-dimensional model of 
adsorption. 

We shall continue this point of view in the present paper, and consider 
the adsorption potential e(r), in the following form 


eo for el 
&(r)-4 €(r) for Pek (4) 
0 for Fe y, 


Above Va, V, and Vy are volumes of adsorbent, of surface phase and of 
free gas phase, respectively. Let us assume later that for reVs, the po- 
tential e has the form 
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where r; is the position on the adsorbent surface, considered as a plane. 
It is like a square-well potential, the y value being simply the slow-collision 
distance between admolecule and adsorbent surface. By choosing the ad- 
sorption potential in the form as given in eq. (5) we assume that he- 
terogeneity does not change the interactions depths y and A, but only the 
strength of the gas—solid interactions. Now, we would like to explain the 
reason, for which such an assumption seems to be a good approximation. 

Let us assume for instance the adsorbate molecules uniformly dis- 
tributed in the adsorbent volume V. Then, in the case of identical ad- 
sorbent molecules interacting via LJ(m—n) potential with an admolecule, 
one gets the following results [1] 


| -3 nis TE * xA n-3 3 m-3 
E oe € [3 - (2 Ji 


where 


E. 2TTN, c" (s* y (7) 
€ mim] — 


In the above x, is the distance from adsorbent surface, n, m, x*, o* and 
s* are parameters in the adsorbent-adsorbate interaction, and r; is the 
appropriate intermolecular distance, N, being the number density (in 
atom/cm?) of the adsorbent, and 
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It is well-known that the functions of LJ(m—n) type are not very sensitive 
to the choice of n and m. The interaction depth in the gas-solid interaction 
potential should therefore depend predominantly upon the value of z*, 
and through it upon s*. Next, it seems that, while the parameters describ- 
ing the attraction strength between the pair of gas—solid molecules may 
change considerably the other parameters describing the slow-collision 
distances, change to a lower degree. For instance, it has been found by 
Crowelland Steele that for rare gases interacting with graphite, 
the interaction between an individual gas-solid pair, can be presented by 
a potential of the form 


w(rz) =-Ą [e j^ Hs Sj. (ra) | (9) 
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From the values tabulated by Crowell and Steele, one sees, that, 
for example, s* value for pair C—CH, is about 10% higher than that for 
pair C—Ne, whereas A value is higher by about 350%. 

The heterogeneity of adsorbents arises from the presence of different 
molecules in the adsorbent volume V;, or from the presence of molecules 
groupped in different ways. An interaction between an admolecule and 
a group of adsorbent :nolecules may be approximately treated as an in- 
teraction with some new chemical individuals. In both cases, the slow- 
-collision parameters should not change considerably, when compared 
with the changes in attractive forces. Since the interaction depths, as 
pointed out, depend predominantly upon the slow-collision parameter 
s*, hence the heterogeneity does not affect remarkably the interaction 
depth A. The heterogeneity effects should appear as rather different attrac- 
tion energies in the surface phase V;. 

Now let us return to the third gas-solid virial coefficient. When eva- 
luating the integral (3), we shall neglect contributions from all other ad- 
molecules configurations, except for the ,,quasi-planar" configurations, in 
which both admolecules are within surface phase Vs, ie. in which both 
0<2,SA and 042,~A. According to the discussion given above, and ac- 
cording to the results of the papers [6, 7], the contributions from the con- 
figurations in which z,. z;€V, are of order 10!? higher than those in which 
admolecules are out of surface phase Vs. 

The theoretical investigations of B4,, for homogeneous surface based 
on two-dimensional models of adsorption, lead to satisfactory agreement 
with the experiments. Consequently the deviations from planarity in 
adsorbed molecules are very small. This is possible when the interaction 
depth A is small. Thus, interaction energy w; can be considered with good 
approximation as a function of intermolecular distance r;; only, between 
the adsorbed molecules [14] 
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It means that the slow-collision diameter remains practically unchanged, 
and therefore we shall take later s=s. 
Next, we shall accept this type of interaction for the case of hetero- 


geneous surfaces, by introducing the concept of an "effective parameter 
ejf”, being a kind of average value of wy averaging over all possible con- 


where: 
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figurations of two adsorbed molecules on the adsorbent surface. With this 


concept we have 
6 
Wz (Mp) ») = kWgp|(7 = TU = (&) (11) 


Let us consider the two spatial distributions: 1) the fully patchwise spatial 
distribution, 2) the fully random spatial distribution. The patchwise spatial 
distribution means, that for the -najority of configurations, in which two 
admolecules are at interaction distances, the both admolecules have the 
same sie energy. In such a case eq. (3) may be written in the form 
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The random spatial distribution means that if one admolecule is located 
anywhere on the adsorbent surface, the second one will have an average 
energy e, equal to 


E - f © X (0)de' (13) 
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where Q is the range of possible variations of e. In this case the eq. (3) 
may be rewritten in the form 


co 


is d ! i j 
mos op(ór | [eli e]ue . 
0 
"fh p 42 6 
. 2a [7 exp - pelt E) E) | -14 dr, (14) 
o 


The two spatial distributions discussed above are some ideal cases, the 
real spatial distribution being only their better or worse approximations. 
The eqs. (12) and (14) provide a method of detection, if an investigated 
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adsorption system obeys the patchwise spatial distribution rather than the 
random one. For this purpose we notice that both the integral and 
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the average e can be easily calculated, by the method developed 
in the first part of this paper. Thus, the expression for B4, from eqs. 
(12) and (14) can be regarded as two-parameter function of tempera- 
ture, with unknown parameters A and ©e;. These parameters may be 
found when fitting best the expressions (12) and (14), with respect 
to experimental values of B;s. Although parameters A and Wey; are chosen 
optimally the expressions (12) and (14) will approximate the experimental 
values of B4, to a various degree. By comparing the experimental de- 
viations one may detect whether the expression (12) or (14) is more ade- 
quate to describe the experimental behaviour of B3, and consequently, 
either the patchwise or the random spatial distribution is a more appro- 
priate model for the investigated system. 


NUMERICAL RESULTS AND DISCUSSION 


In order to demonstrate the applicability of the procedure proposed 
above, we shall consider the experimental results of Sams et al. [15] 
concerning the adsorption of argon on graphitized carbon P33 (2700). 

As it was mentioned in the first part, these are most accurate experi- 
mental data reported thus so far. It will be interesting to compare the 
results obtained here with the earlier theoretical results [12, 13, 16, 17], 
obtained with the assumption that the graphite P33 (2700) is a homoge- 
neous adsorbent. Table 1 presents the results of our calculations, per- 
formed on the computer ODRA-1204 in the UMCS Department of Nu- 
merical Calculations. The numerical results show that the assumption of 
the patchwise spatial distribution of adsorption energies leads to negative 
values for both A and wej which are of course physically meaningless. 
The physical meaning of the parameters A and wes; suggests that they 
must have still positive values. 

Now, let us consider the values for A and wesp, obtained by assuming 
the fully random spatial distribution. From the values tabulated by 
Hirschfelder et at. [18] one sees, that parameter w for two argon 
atoms in bulk phase interacting via LJ (12—6) potential is about 240 cal/ 
(mole. Johnson and Klein [17] assumed that the surface of graphite 
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Table 1. The results of numerical calculations for the experimental results of Sams, 
Constabaris and Halsey 


The assumed The gas-solid The gas-solid 
model interaction interaction 
of adsorption depth A energy ert 
Patchwise d = V 
distribution 2.65 10-5 A 478 cal/mole 
Random 
distribution 1.60 A 226 cal/mole 


P33 (2700) was fully uniform, and accepted the Sinanoglu-Pitzer potential. 
They obtained the value w= 220 cal/mole and s=3.45 A, whereas in bulk 
phase 3.41 A. The theoretical analysis of Wolfe and Sams [12] yields 
the values 220 cal/mole and 3.43 À, respectively, when accepting the Si- 
nanoglu-Pitzer interaction potential, and 118.8 cal/mole and 3.46 À, when 
accepting the LJ (12—16) potential. The small differences, which appear 
are, in our opinion, due to heterogeneity of the surface. The above results 
show that the spatial energy distribution for argon on graphite P33 (2700) 
is a highly random one. The best-fit results for the effective slow-collision 
diameter s, obtained by Johnson and Klein, and Wolfe and Sams, justify 
its identification with appropriate value for bulk phase, made in our cal- 
culations. Both Johnson and Klein, and Wolfe and Sams have found that 
the calculation of surface area for this adsorption system is insensitive 
to the three-body effects. According to theoretical results of Steele [3] 
and Hill [4], such a situation is characteristic of the adsorption systems 
with random energy distribution. Now let us consider the interaction 
depth A. It is equal to 0.47s— 1.60 A. Thus, it is in fact a small parameter, 
as it was assumed by us. 


On the basis of the above discussion we may conclude that the distri- 
bution of adsorption centers on graphite P33 (2700) surface must be de- 
scribed by the random distribution function. 
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STRESZCZENIE 


W części pierwszej pracy wskazano metodę wyznaczania funkcji dystrybucji 
energii adsorpcji z temperaturowej zależności drugiego współczynnika wirialnego 
gaz—ciało stałe. W niniejszej części wskazano, w jaki sposób można z temperatu- 
rowej zależności trzeciego współczynnika wirialnego gaz—ciało stałe uzyskać do- 
datkowe informacje na temat przestrzennego rozkładu energii adsorpcji. 


PESIOME 


B neppow uacTru paOoTbI Obi] yKa3aH HOBbIM MeTOĄ orpegezenusa uWdcdepeu- 
IHMaJIbHOrO pacrpeąejleHua SHepruu agcopõunnu M3 BTOporo BMpMaJlbHoro Kosddwu- 
nMeHTa ra3—TBepyoe TeJo. B cjenyrolień uacTM IIOKA3AHO, KAKUM OÓpa30M MO2JKHO 
HOJXy4MTb ZaJbHeAnnune unHdopMaHu4M Ha TeMy IpOCTpaHCTBeHHOTO pacrpejleJeuus 
3HepTMM AĄCOPÓNMKM 43 TpeTbero BUPMAJIbHOTO Ko»dduuueura. 
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On. the basis of the above discusion We may conclude th wt the distri- 
tution of adaorptíon centers on grzpblig P38 | 1200) surface „państ bu des 
scribed by the random distribution funy 
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